Spatial Persistence of Fluctuating Interfaces 
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We show that the probabihty, Pail), that the height of a fluctuating (d+ l)-dimensional interface 
in its steady state stays above its initial value up to a distance I, along any linear cut in the d- 
dimensional space, decays as Po(0 ~ l~ ■ Here is a 'spatial' persistence exponent, and takes 
different values, 6s or ^o, depending on how the point from which I is measured is specified. While 
6s is related to fractional Brownian motion, and can be determined exactly, 6o is non-trivial even 
for Gaussian interfaces. 



PACS numbers: 05.70.Ln, 05.40. +j, 02.50.-r, 81.10.Aj 



Fluctuating interfaces have played the role of a 
paradigm in nonequilibrium statistical physics for the 
last two decades, with applications ranging from growth 
problems to fluid flow [y . Traditionally the different uni- 
versality classes of the interfaces are specified by the ex- 
ponents associated with the dynamic correlation func- 
tions. Recently the concept of persistence |^, i.e. the 
statistics of first-passage events, was used to characterize 
the temporal history of an evolving interface p] . Starting 
from an initially flat interface, it was found that the prob- 
ability that the interface at a given point in space stays 
above its initial height up to time t decays as a power-law 
in time with an exponent 9}^™ that is nontrivial even for 
simple 'Gaussian' interfaces (i.e. with dynamics governed 
by a linear Langevin equation) |g,Q| . On the other hand, 
the statistics of returns to an initial profile chosen from 
the ensemble of equilibrium configurations was shown to 
be governed by a different exponent 0*<^™ . The exponent 
Qtem .^^g shown be the identical to the first passage ex- 
ponent of the fractional Brownian motion (fbm) which is 
known exactly |g|]. In |3| these two regimes were termed 
the coarsening and stationary regimes respectively. 



The exponents O}^™' and 6'*'^™ describe first-passage 
properties in time. A question that naturally arises is: 
what about the first-passage properties in space? Is there 
a power-law distribution and associated nontrivial expo- 
nents for the spatial persistence of an interface? In this 
paper we address this important question and show that 
indeed in the steady state (long time limit), the proba- 
bility Po(0 that an interface stays above its initial value 
over a distance I from a given point in space decays as 
Pq{1) ^ l^^ for large I. Furthermore, analogous to its 
temporal counterpart, the spatial persistence exponent 6 
again takes different values, 9q and 9 a, according to the 
"initial" condition, in this case specified by the height 
and its spatial derivatives at the point (the "initial point" 
in space) from which I is measured. If the initial point 
is sampled uniformly from the ensemble of steady state 
configurations, the relevant exponent, 9s, is related to the 
Hurst exponent of the spatial roughening. Conversely, if 



the initial point is such that the height and its spatial 
derivatives are finite, independent of the system size (we 
will call this a finite initial condition), the corresponding 
exponent ^o is a new, nontrivial exponent. 

A second interesting question asks: Is there a morpho- 
logical transition in the stationary profile of a fluctuating 
interface if one changes the mechanism of fluctuations by 
tuning either the dynamical exponent z or the spatial 
dimension dl This question has important experimen- 
tal significance. For example, recent Scanning Tunneling 
Microscope (STM) measurements have shown [p| that the 
dynamical exponent z characterizing the fluctuations of 
single layer Cu(lll) surface changes from z = 2 at high 
temperatures to z = 4 at low temperatures ^^. Thus 
by tuning the temperature, and hence z, the STM mea- 
surements may possibly detect a morphological transition 
if there is one. We show here that Gaussian interfaces in- 
deed exhibit a morphological transition across a critical 
line Zc{d) = rf + 2 in the {z — d) plane as one tunes ei- 
ther z or d. For z > Zc{d), the steady state profile is 
smooth with a finite density of zero crossings along any 
linear cut in the d-dimensional space. On the contrary, 
for z < Zc{d), the density of zero crossings is infinite and 
the locations of zeros are nonuniform and form a fractal 
set. 

A seemingly unrelated problem is the usual temporal 
persistence of the stochastic process B , 






(1) 



where x is the position of a particle and ?y(i) is a Gaus- 
sian white noise with zero mean and delta correlation, 
{v{t)ri{t')) = S{t - t'). The probability, Po{t), that x 
does not change sign up to time t decays as Po{t) ~ t"^*^") 
for large t, where the persistence exponent 9[n) depends 
continuously on n H . The exact value of the exponent 
is known only for n = \ (the usual Brownian motion), 
9{1) = 1/2 iiq], and n ~ 2 (the random acceleration 
problem), 61(2) == 1/4 ||ll|,|l§. The latter problem has re- 
cently attracted attention in connection with the problem 
of inelastic collapse in granular materials |]l3[ . 



In this paper, we establish an intricate mapping be- 
tween the process in Eq. (P and the steady state profile 
of Gaussian interfaces evolving via the Langevin equa- 
tion, 



dh 
'dt 



^^ = -i-y'f'h + C 



(2) 



where V^ is the d-dimensional Laplacian operator, z is 
the dynamical exponent associated with the interface and 
^(r, t) is a Gaussian white noise with zero mean and 
{^{f,t)^ir',t')) = 26{r~r')5{t-t'). The continuum equa- 
tion (g) is well defined only for z > d. For z < d, one 
needs a nonzero lattice constant in real space or an ultra- 
violet cut-off in momentum space. Below we will assume 
z > d. 

Let us summarize our main results. We find two inde- 
pendent spatial persistence exponents 6q and 6s depend- 
ing on the initial point in space (from where measure- 
ments start). The two main results are: (i) By exploiting 
the mapping to the stochastic process (0), we show that 
for a 'finite' initial starting point. 



9{n), where n= {z - d+ l)/2. 



(3) 



and (ii) For initial points sampled from the steady state 
configurations, however, we show exactly that 



(3/2-n, l/2<7i<3/2 
\0, n>3/2. 



(4) 



where n = {z ~ d+ l)/2. 

Exploiting the two exact results, 6{1) — 1/2 and 
e{2) = 1/4 in Eq. (|), we find that for any {d + 1)- 
dimensional interface with z — d = 1, 9q = 1/2 and for 
z — d = 3, 9q = 1/4: exactly. An explicit example of the 
former case is the (1 + l)-dimcnsional EW model with 
z = 2 and d = I, giving ^o — 1/2 Q. An example of the 
latter case is the (l-l-l)-dimensional continuum version of 
the Das-Sarma Tamborenea (DT) model ||l^ with z — 4 
and d = 1, which gives 9o = 1/4, a new exact result. 
For general n > 3/2, we show that the exponent 9{n) 
and hence 9o can be estimated rather accurately by us- 
ing an 'independent interval approximation' (IIA). Thus 
Eq. (H), besides giving accurate, and in some cases ex- 
act, results for the new exponent 9o, also give a physical 
meaning to the stochastic process in Eq. (|l|) for general 
n. We will use this correspondence to further establish a 
morphological transition for Gaussian interfaces, across 
the critical line Zc{d) — d + 2, as one tunes z or d in the 
{z — d) plane. 

Let us consider the n-th derivative of the scalar field 
h{r,t) with respect to any particular direction, say xi, 
in the d-dimensional space, g{r,t) = d"h/dxi. From 
Eq. (H), it follows that the Fourier transform, g{k,t) = 



l9{r,t), 



ik.r^d 



d r then evolves as, 



dt 



^\k\^gik,t)+^,{k,t) 



(5) 



where |fc|^ = fc^ -I- fc| -f . . . and ^i is a Gaussian 
noise with zero mean and correlator, (^i(fc, t)^i{k' , t')) — 
2\ki\'^''5{k+k')5{t-t'). Using Eq. (|) , one can easily eval- 
uate the correlator, {g{k,t)g(k' ,t)) — G{k,t)S{k + k'), to 
obtain, in the steady state, G{k) — |fcip"/|A:|^. Invert- 
ing the Fourier transform gives the stationary, real-space 
correlator, 



{g{xi,X2,---)g{x'i,X2,---)) 



Ifcil 



\k\^ 



Jki(xi-x'i)^d^ 



(6) 



It is now easy to see from Eq. ^ that with the choice 
n ~ {z~d+ l)/2 the integral yields a delta function, i.e., 

{g{xi,X2,X3,...)g{x'i,X2,X3,...)) = DS{xi -x'l) (7) 

where D is just a dimension-dependent number. 

Since the basic Langevin equation, Eq. (H), is a linear 
equation and the noise ^ is Gaussian, clearly /i is a Gaus- 
sian process and therefore its n-th derivative, g, is also 
a Gaussian process. A Gaussian process is completely 
specified by its two-point correlator. Hence in the steady 
state limit i — > oo, one can write the following effective 
equation, which gives the xi-dependence of g at fixed 

X2,X3,. .. 



g{xi,X2,..-) 






'7(a;i), 



(8) 



where ry is Gaussian white noise with zero mean and 
{ri{xi)r]{x'i)) = 5{xi — x'l). Note that we have rescaled 
the xi axis to absorb the constant D. This effective 
stochastic equation will generate the correct two-point 
correlator as given by Eq. (|7|) and also all higher order 
correlations as well since 5 is a Gaussian process. Making 
the identification, xi <^=> t and h <=^ x, Eq. (H) immedi- 
ately reduces to the process in Eq. (mj . Hence the steady 
state profile of a Gaussian interface evolving via Eq. (0) 
can be effectively described by the stochastic equation 
(m with n = {z — d + l)/2. Therefore, the spatial per- 
sistence of the Gaussian interface also gets mapped onto 
the temporal persistence of the stochastic process (Ul). 

Consider the stochastic process (|l|) starting from the 
initial condition at t = where x and all its time deriva- 
tives up to order (n — 1) vanish. It turns out that the 
statistics of first passage events of this process depends 
crucially on whether one starts measuring these events 
right from t = or or if one first waits for an infinite time 
and then starts measuring the events. This is similar to 
the 'coarsening' versus 'stationary' regimes as studied in 
[p|. This can be quantified more precisely in terms of 
two-time correlations. Integrating Eq. (1) n times gives 



x(i)- 



1 



r(n) J, 



r^it,)it-t,r-^dh, 



(9) 



where T{n) is the usual Gamma fmiction. From Eq. (|9|), 
which provides a natural continuation from integer n to 
real n, it is easy to evaluate the correlation function, 

Cit, t', to) - {[x{to + t) - x{to)] [x{to + t') - x(to)]) 
= A{to + i, to + t') - A{to + 1, to) 
-.l(to,to + t') + A(to,to), (10) 

where the autocorrelation A{t^f) is given by, 

1 p'min{t,t') 



A{t,t') 



r'in) Jo 



(t-ti)"-i(t'-ti)"-Mti. 



(11) 



Thus the correlation function that fully characterizes the 
Gaussian process and thereby its persistence, depends 
explicitly on the 'waiting' time toi the starting point of 
measurements. It turns out, however, that there are only 
two asymptotic behaviors controlled respectively by the 
to = and tg —f oo fixed points. Any finite to flows into 
the to = fixed point. 

To see how the two exponents, Oq and 6s-, emerge, con- 
sider first the limit to — > oo. From Eq. (|l^), we get after 
some algebra, 



the exponent Oq is nontrivial. To find it one is con- 
fronted with the following problem: given a Gaussian 
stationary process X{T) with a prescribed correlator 
C{T) = {X{0)X{T)), what is the probability Po{T) that 
the process X{T) does not cross zero up to time T? 

For general correlator C{T), Po{T) is hard to solve 
]l8|. For n = I, we have, from Eq. (|l|), Ci{T) = 
exp(— r/2), a pure exponential. This corresponds to a 
Markov process, for which the exact result is known: 
Pq{T) - exp(-T/2) - t-i/2 implying 6lo(l) = 1/2. For 
higher values of n, the process X{T) is non-Markovian 
and hence to determine 6'o(n-) is harder. Fortunately for 
n = 2 where C2{T) = | exp(-T/2) - i cxp(-3T/2), the 
exact result is also known: Po{T) ~ exp(— r/4) ~ t~^^'^ 
|ll[. For general n, one can extract some useful informa- 
tion about the stochastic process by studying the short- 
time properties of the correlator C„ (T) . Expanding Eq. 
(|l|) for small T we get, 



C(t,t',oo)^a„[t2"-i+t'2"-i-|t-tf"-i], l<n< ^, 



3 

2' 

(12) 



where a„ = -r[n]r[l - 2n]/r[l - n] > and /3(n) = 
[n^ - n + 1 - (n - 2)/(2r7. - 3)]/(2n - 1) > 0. To calculate 
the persistence properties from this correlator, we note 



that the form of the correlator for 



<n< 



in Eq. 



^) 



is precisely that of fbm with Hurst exponent H = n~ ^. 
The corresponding persistence exponent is known exactly 
§,|6| to he Bs = 1 - H = |-n= {d + 2 - z)/2, 
from Eq. (||), provided z < d + 2. For z > d + 2, 
i.e. n > 3/2, Eq. ([l2[) gives the normalized correlator 
C{t, t', cx))/VC(t,t,oo)C(t',t',cxj) = 1, implying 0s = 0. 
We now turn to the limit to = 0. From Eq. (10), 



we find that the process is non-stationary in time since 
the correlator depends on both t and t' and not just 
on the difference \t — t'|. However, using the standard 
transformation ig], X = x{t) / yJJxHff) and T = Int, 
one finds that X{T) becomes a Gaussian stationary pro- 
cess in the "new" time variable T, with a correlator 
C„(T) = {X{Q)X{T)), 

Cn{T) = (2 - i) e~^/2 p^^ _ „^ 1- 1 + „. g-T)^ (^3) 
n 

where F{a, b; c; z) is the standard hypergeometric func- 
tion [O. This form of the correlator suggests that 



Cn{T) 




1/2 <n< 3/2, 
n = 3/2, 
n > 3/2, 



(14) 



where a„ = r(n)r(2 - 2n)/r(l - n). Thus for n > 3/2, 
the process is "smooth" with a finite density of zero cross- 
ings that can be derived using Rice's formula UM, 



v^cm 



2n-l 



1 



27r V 2n - 3 ' 



(15) 



where C„"(0) is the second derivative at the origin. For 
1/2 < n < 3/2, the density is infinite and the zeros 
are not uniformly distributed but instead form a frac- 
tal set with fractal dimension dj = n — 1/2 ||l^. The 
case n — 3/2 is marginal, and the integral in Eq. ( |l^ ) 
can be evaluated in closed form, C3/2(T) = cosh(r/2) -t- 
sinh^(T/2) ln[tanh(T/4)], with the density still divergent 
but only logarithmically. A physical example that corre- 
sponds to this marginal case n = 3/2 is the 2-d DT model 
with z — A. The exponent 0{n) diverges as n — > 1/2 from 
above. For n < 1/2 or cquivalently z > d the continuum 
equation is not well defined and one needs a non-zero lat- 
tice constant in space (for the interface problem) or time 
(for the stochastic process in Eq. (P). 

This transition at n — Uc = 3/2 can now be di- 
rectly interpreted in the interface context using the re- 
lation (0). Gaussian interfaces with finite initial condi- 
tion undergoes a morphological transition at Zc — d + 2. 
For z > Zc{d) {n > 3/2), the surface has a finite den- 
sity of zero crossings, given by (0) and dlq) as p = 
(l/27r)-y/(z — d)/{z — d — 2), along any linear cut in the 
d-dimensional space. On the other hand, for z < Zc{d), 
the surface crosses zero in a nonuniform way and the den- 
sity is strictly infinite. If it crosses zero once, it crosses 
subsequently many times before making a long excursion. 



This transition is reminiscent of the "wrinkle" tran- 
sition found recently by Toroczkai et. al. [g] in (1 + 1)- 
dimensional Gaussian interfaces. However there is an im- 
portant difference. Toroczkai et. al. studied the density 
of extrema of a Gaussian (1 -|- l)-dimensional surface and 
found that the density of extrema is finite for z > 5 and 
infinite for z < 5. Their results can be understood very 
simply within our general framework by noting that an 
extremum of the surface h corresponds to dxh = 0. Thus 
the density of extrema of the surface h corresponds to the 
density of zeros of the derivative process dxh. One can 
then follow our chain of arguments mapping a Gaussian 
interface to the stochastic process in Eq. (|l|) and one finds 
that the derivative process also maps to Eq. (|l|) , but with 
n = {z — d — l)/2, i.e. n gets replaced by (n -I- 1) in Eq. 
(g|) due to the extra derivative. Using ric = 3/2 again, we 
find that the transition for the derivative process occurs 
a,i z ~ Zc = d + A. This is thus a generalization of the 
d= 1 result (zc = 5) of Ref. §. 

For n > 3/2, where the process is smooth, one can 
apply the IIA which assumes that successive intervals 
between zero-crossings of the process X(T) are statis- 
tically independent. The IIA was used very successfully 
for rather accurate analytical estimates of the persistence 
exponent for the diffusion equation ||9|,po||. According to 
this approximation, the exponent 9 for a general smooth 
Gaussian process with correlator C(T) is given by the 
first positive root of the following transcendental equa- 
tion i, 



d < z < d+2, and C{x) ~ x'^L^~'^~'^ for z > d+2, where 
the system size L acts as a regulator for the small-fc di- 
vergence {kjnin ^ 1/-^) in the latter case. This analysis 
also illustrates the difference between finite and steady 
state "initial" points: if h{0) is fixed, i.e. not averaged 
over the steady state distribution, one obtains, for large 
\x\, C{x) ^ |a;|^~'^, independent of L, for all z > d. The 
same result follows from Eq. ffl) after the replacements 
X ^ h, t ^ X. 

In summary, we have shown that fiuctuating interfaces 
exhibit a form of spatial persistence analogous to the tem- 
poral persistence exhibited by stochastic processes. For 
Gaussian interfaces the analogy is precise - the spatial 
fluctuations in the stationary state are isomorphic to the 
stochastic process (nl), with n given by Eq. (0). A non- 
Gaussian process for which exact results are possible is 
the (1 + l)-dimensional KPZ equation, since the station- 
ary probability distribution of the interface fleld is given 
by, P{{h{x)}) - exp[- J{dxh)'^dx] §. Thus the station- 
ary interface can be described by the effective Langevin 
equation, dxh = vi^): the n = 1 version of Eq. (|l|). From 
the exact results 6(1) — 1/2, and H = 1/2 we immedi- 
ately obtain 9q — 9s = 1/2, in agreement with Ref. [[l4[ . 
For the KPZ equation in higher dimensions, and other 
non-Gaussian interface models, the determination of 0Oi 
in particular |22|] , remains a challenge. 

We thank M. Barma, G. Manoj, M. Marsih and C. 
Castellano for useful discussions. 
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where p = y/—C"{0)/TT is the density of zero crossings. 
Applying this formula to our problem, with Cn{T) given 
by Eq. (13) and p by Eq. dlq), we can obtain estimates 
for 9{n) for arbitrary n > 3/2. 

For example, putting n = 2 in Eqs. (O) and (16), we 
get 9{2) — 0.26466 . . . (IIA) which can be compared to 
the exact result 9{2) — 1/4. Similarly for n = 3, we 
flnd 6'(3) = 0.22283 . . .(IIA), to be compared to the value 
9(3) « 0.231 ± .01 obtained by numerical simulation pi[ . 
The limit n — > cxd is rather interesting. From Eq. (13), 
we get, Coo{T) = scch(T/2) which happens to be also 
the correlation function for the two-dimensional diffusion 
equation P,pO|. For this, one can also obtain an IIA esti- 
mate, 9{oo) = 0.1862 . . . (IIA) [||,|0| which is again rather 
close to the numerical value 9{oo) k, 0.1875 ± 0.0010 [|. 
Thus for all n > 3/2, the analytical IIA estimates for 
9{n) seem to be rather accurate. 

In the context of a stationary interface, the morpho- 
logical transition at n = 3/2 (i.e. z — d + 2) \s associ- 
ated with the familiar transition from a rough to a super- 
rough interface, extracted from the mean-square height 
difference C{x) = {[h{x) — /i(0)]^). Using the station- 
ary correlator, (/ik/i-k) ^ l/|k|^ gives C{x) ~ |a;|^~'^ for 



[4; 

[6; 
[7: 
[s; 
[9; 
[lo; 

[11 



For a recent review, see J. Krug, Adv. Phys. 46, 139 

(1997). 

For a recent review on persistence, see S.N. M ajum- 



dar, Curr. Sc i. 77, 370 (1999), also available on cond 



mat/9907407 



J. Krug, H. Kallabis, S.N. Majumdar, S.J. Cornell, A.J. 
Bray, and C. Sire, Phys. Rev. E 56, 2702 (1997). 
H. Kallabis and J. Krug, Europhys. Lett. 45, 20 (1999). 
M. Giesen and G.S. Icking-Konert, Surf. Sci. 412/413, 
645 (1998). 

Z. Toroczkai and E.D. Williams, Phys. Today 52, 24 
(1999). 

S.V. Khare and T.L. Einstein, Phys. Rev. B 57, 4782 
(1998). 

Z. Toroczkai, G. Korniss, S. Das Sarma and R.K.P. Zia, 
Phys. Rev. E 62, 276 (2000). 

S.N. Majumdar, C. Sire, A.J. Bray, and S.J. Cornell, 
Phys. Rev. Lett. 77, 2867 (1996). 

W. Feller, Introduction to Probability Theory and its Ap- 
plications, 3rd ed. Vol-1 (Wiley, New York, 1968). 
T.W. Burkhardt, J. Phys. A 26, L1157 (1993); Y.G. 
Sinai, Theor. Math. Phys. 90, 219 (1992). For a par- 
ticularly simple derivation, see also T.W. Burkhardt, 
preprint (2000). 



[121 G. De Smedt, C. Godreche, and J.M. Luck, cond 



mat/0009001 



[13] 



[14] 
[15] 

[16] 



S.J. Cornell, M.R. Swift, and A.J. Bray, Phys. Rev. Lett. 

81, 1142 (1998); M.R. Swift and A.J. Bray, Phys. Rev. 

E. 59, R4721 (1999); T.W. Burkhardt, J. Franklin, and 

R.R. Gawronski, Phys. Rev. E. 61, 2376 (2000). 

D. Das and M. Barma, Phys. Rev. Lett. 85, 1602 (2000). 

S. Das Sarma and P. Tamborenea, Phys. Rev. Lett. 66, 

325 (1991). 

A. Hansen, T. Engoy, and K.J. Maloy, Fractals 2, 527 

(1994); M. Ding and W. Yang, Phys. Rev. E 52, 207 

(1995). 



[17] LS. Gradshteyn and I.M. Ryzhik, Table of Integrals, Se- 
ries, and Products, p-1040 (Academic, Orlando, 1980). 

[18] I.E. Blake, and W.C. Lindsay, IEEE Trans. Inf. Th. 19, 
295 (1973). 

[19] S.O. Rice, Bell. Syst. Tech. J. 23, 282 (1944). 

[20] B. Derrida, V. Hakim, and R. Zeitak, Phys. Rev. Lett. 
77, 2871 (1996). 

[21] S.J. Cornell, private communication. 

[22] For (2 + l)-dimensional KPZ interface, the relation 8s = 
1 ~ H has recently been verified by G. Manoj, private 
communication . 



